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1 Introduction 

Since the observation by Glover-Lazer-McKenna IS that a simple harmonic oscillator with a piecewise linear 
stiffness (jumping nonlinearity) contributes to the explanation of 

the failure of the Takoma bridge, the studying of periodic oscillations in such models got a lot of attention 
of mathematicians; see the recent survey HI and the papers |[T2l l2l. Also, new engineering studies of impact 
oscillators open up a large potential for challenging extensions of these results. In fact, Ivanov H argued 
that harmonic oscillators with a jumping nonlinearity with one part of the force field nearly infinite is a better 
model for describing the bouncing ball, rather then its limit version for an impact oscillator. In our modeling 
the resulting system of differential equations is singularly perturbed, but as we discuss below, the classical 
singular perturbation theory does not apply. In this paper we develop an averaging-like approach which solves 
the problem in a weakly nonlinear case. For a discussion of the use of averaging method for regular impacting 
systems we refer to HI. To be explicit, our approach concerns the existence and stability of periodic oscillations 
of the following system 

X + X = ef{t,x,x,e), x > 0, 

where f,g G C^(M x M x M x [0, l],]R),e > is a small parameter, ^ G M as e ^ 0. System [TT] 
can be considered as a smoothed version of a system with impacts. We will study resonance oscillations and 
assume, therefore, that 

f{t + vr, u, V, e) = f{t, u, v, e),g{t + vr, u, v, e) = g{t, u, v, e). 

System (11.11 ) represents a natural singular perturbation description of impact phenomena which is different 
from the usual approaches. Our main result (Theorem 1) states that the emergence of asymptotically stable 
TT-periodic solutions in (11.11) from vr -periodic cycles of non-smooth limiting the system 

X + X = 0, X > 0, 

x{t-0) = -x{t + 0), x{t) = 0, 

can be studied by a special form of the averaging method combined with a suitable scaling of time when 
solutions pass the half plane x < 0. This involves the use of the implicit function theorem for a non-smooth 
problem in the limit as e — > 0; this is possible by introducing a suitable Poincare map. The result is a change 
of stability of a fixed point when its eigenvalues enter the unit disc from outside through the imaginary axis. 



2 Main result 

We prove the following theorem. 
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Fig. 1: Illustration of the change of variables i2.5t -l l277l l. 



Theorem 2.1 Let f,g£ C'^{R X M X M, M) be Tr-periodic with respect to time. Define 

iin(r + 6 
cos(r + 



_ f-V sin(r + 0) \ 

P{A,e) = - / 1 \{f{T,Acos{T + 6),-As{n{T + d),Q) -2loqAcos{t + d))dT- 



\^ 
TT / sin(T + 6' + 7r) \ 
-/ I , , (/(r,ylcos(T + e + ^),-^sin(T + + 7^),O)- 

-2ojqA cos(t + 6 + ■K))dT - ujq] i^^^ ^ 2 ) ^ 5 (^^ - 61, 0, sin (^s + , 0^ ds. 

If P{Aq, 9q) = Ofor some {Aq, 9q) ^^M. x (0, vr) and the real parts of eigenvalues of p'{Aq, Oq) are negative, 
then, for any e > sufficiently small, equation di.il ) has a unique n-periodic solution satisfying 

^eit) xo{t) as e ^ pointwise on [0,7r] \ |— — 6q^ 
where xq is the unique n -periodic solution of the equation di.2D with initial condition 

(x(0),i;(0)) = (^0 cos 00,-^0 sin 0o)- 
Moreover, the solution is asymptotically stable. 
Proof. Rewrite system (11.11) as follows (see Fig. Ill) 

(1 - euJeY (1 - ewj^ (1 - SLOeY 

so, that any solution of the reduced system (e = 0) 

1 

^ + T, 7^2; = 0, X > 0, 

(1 - euJeV 

1 

(e)2(a;^.)2 

is TT-periodic. Let us introduce variables {A, 9) as follows 

1 



Acqs 



1 - euJe 



A 1 • ^ 1 

-A sm 



1 — ELOp V 1 — 



0,|(l-ea;. 



(2.5) 



2 



eA cos 



-A-^ sin —(0 - -(1 - eu,)) + - 



1 — euj, 



-(l-ea;e),-(l + ea;e 



A cos 



-A- 



1 - 



^(l + eo;,)) +|+7r) , 



1 



1 — eojp 



■ sm 



1 



1 — eujp 



which transforms equations (I2.3l )- (l2.4b to the following system 
'A\ /O 



+ £Gi{t,A,e,£), if 
+ eG2(t,A^,e), if 
+ eG3(t,^,e,e), if 0G 



vr 



0,-(l-ea;, 



-(l-ea;e),-(l + ea;£) 



-(l + ewe),7r 



where 



G'i(t,Ae,e) 



-(1 — euje) sin 



■ COS 



G2{t,A,e,e) 




— (1 - euje) sm 



■9 ( 

-A 



COS 

1 — eujp \ EUJ, 



1 



1 — euJt_ 

■ — (1 -eu)e)^e COS I 

A \ eujp 



5 t, 



1 - 

, 1 

A sm 

1 — ELOp 



COS +^ 

1 /„ vr , ,\ 7r\ 



G3{t,A,9,e) 



-(1 — euj^) sin 
/ ^t, A cos 
-A- ' 



1 



1 - £LJe 
1 

1 — euj^ 



(l+euJe)) + - + vr 



2(1 + ^^^ 



" 1 — ew, 



■ sm 



e 

We 



1 - eto^ 



TT 



IT 

H h vr 

2 

vr 



(l + ewe) +- + vr ,e 



(1 - ELVe 



(2 + ews) cos 



1 



1 — euj^ 



TT 



-{l + eus)) +- + 7r 



TT 



COS 



A 

f ( t,A cos 
1 



1 



1 — eujp 



1 - etor. 



TT 



vr 



-A- 



1 — ew^ 



■ sm 



1 — euJr 



TT 

V-+TT 
TT 

H h vr 

2 

vr 



v 



We 



r(2 + eWe) COS 



1 — eujp 



TT 



TT 



(1 + eUe) + n + ^ 



lit ^ {A{t),9{t) — t) is an asymptotically stable vr-periodic solution of (I2.8l) - (l2.10b then (x,x) defined by 
(I2.5I) - (I2.7I) is an asymptotically stable vr-periodic solution of (I2.3I) - (I2.4I) . To prove the existence of asymptoti- 
cally stable vr-periodic solutions of equations (I2.8I) - (I2.10I) we show that each solution 
(A{-,A,e,e),d{-,A,d,e)) of (irst- dTTOl) with initial condition (A{0,A,9,e),e{0,A,e,e)) = {A,e) is de- 
fined on [0, vr] whenever {A, 9) belongs to a small neighborhood of {Aq, 9q), and that the map 



Pe{A, 9) = {A{tt, a, 9, e),9{TT, A, 9, e) - vr) 



(2.11) 



contracts in this neighborhood. 

Step 1, First we show that solution t 1-^ (A{-, A, 9, e),9{-,A, 9, e)) of (IT8])- (I2.10I ) on [0, vr] can be conse- 
quently sewed by solutions of systems (12.81 ). (12.91 ) and (12.101 ). 

Denote by t 1-^ (Ai{-,to,A,9,e),9i{-,to,A,9,e)), i = 1,2,3, the solutions of (ITSl) . ( |2!9l ). (12.101 ) respec- 
tively with initial condition {Ai{to, Iq, A, 9, e),9i{to, to. A, 9, e)) = (^4, 9). Put 



F,{T,A,9,e) 



1 — eto. 



^ 9,{T,0,A,9,e)-^. 



Since Fi 



TT 



^Oi Aq, 9o,0] — 9i 



TT 



9o,0,Ao,9o,0 



T 



TT 



TT 



TT 



(1) V2 



and 



),0,Ao,^o,0 =1 



then by the implicit function theorem E Ch. X, § 2, Theorems 1 and 2] there exists Ti € C^(E x M x [0, 1], M) 
such that 

Fi{Ti{A,9,e),A,9,e)=0, \A-Ao\<6, \9 - 9o\ < 6, [0,<5), 
where S > sufficiently small. Or, equivalently, 



1 — eto, 



^ 9i{Ti{A,9,e),0,A,9,e) 



2' 



\A-Ao\ < 6, \9-9o\ <6, e£ [0,5). 



Therefore, the solution of system (12.81 ) with initial condition {A, 9) at t = approaches the threshold of 
switching to ( |2!9l ) at time Ti{A, 9, e). 
Now we show that the solution 



vr 



A2 [■,Ti{A,9,e),Ai{Ti{A,9,e),A,9,e),-{l-eu:e),e) , 
92 (•, Ti {A, 9, e),Ai (Ti (A, 9, e),A, 9,e),^(l- euj,),e 



stays till some time T2{A, 9, e) in 



[0,oo)x ^{l-euJe),^{l+euJe 
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and that T2{A, 6, e) is given by 

T2{A,e,e) = Ti{A,e,e)+ef2(A,9,e), 
where f2 G Ci(M x R x [0, 1], M). To do this consider 

F2{T,A,6,e) = —{e^{Ti{A,9,e)+eT,n{A,6,e), 



Ai{Ti{A, 9, e),0, A, 9, e),9i{Ti{A, 9, e), 0, A, 9, e),e) - |(1 - eo;,)) - vr. 



F2{T,A,9,e) = { 



^ — {O2 (Ti {A, 9,e) + eT, Ti {A, 9,e),Ai (Ti {A,9,e),A,9,e), 
^{1 - £uJs),e) -^{1 - euJe)^ - TT, if e > 0, 
— T - vr, if e = 0. 

UJQ 



Let us convince ourself that the function F2 verifies the assumptions of the implicit function theorem ||5l Ch. X, 
§ 2, Theorems 1 and 2] at the point (T, A, 9, e) = (wott, ^Oj Since 

|(1 - EiVe) = 92 (r.iA, 9, e),T,{A, 9, e),MTi{A, 9, e),A, 9, e), |(1 - eo;,), 

then 

limF2{T,A,9,e) = lim —(92yn){Ti{A,9,e) + X{A,9,e)eT,Ti{A,9,e), 

MTi{A, 9, £),A, 9, £), ^(1 - eujs),£) eT - vr = — T - vr, 
that is F2 is continuous at e = 0. Here X{A, 9, e) € [0, 1]. Furthermore, we have 

(^2)t - ^0 + vr, Ao, 9o, o) = — / 0. 

Therefore, the implicit function theorem lH Ch. X, § 2, Theorems 1 and 2] allows as to conclude that there 
exists f2 G Ci(IR X M X [0, 1], M) such that 

— (MTiiA, 9, e) + ef2{A, 9, e), Ti{A, 9, e), Ai{Ti{A, 9, e), A, 9, e), 
^il-eLOe),e)-^il-eLOe)) =tt, \A-Ao\<6, \9 - 9\ < 5, eG [0,5). 



2' " ' 2 

Since ^0 £ (0, vr), then 5 > can be diminished in such a way that 

^3(vr, ri(A, 9, e) + ef2{A, 9, e),A2{Ti{A, 9, e) + ef2{A, 9, e), 

Ti{A,9,e),Ai{TiiA,9,e),A,9,e),^{l-eu;e),e),^{l+euJs),e) < 
vr 

< — — 9 + TT + Tr, fore>0 sufficiently small 

and any |j4 — Ao| < 6,\9 — 9o\ < 6. Therefore, the solution of (12.101 ) with initial conditions under consideration 
does not meet the line of discontinuity during time {Ti{A, 9, e) + eT2{A, 9, e), vr]. 



5 



. {A2iTi{A,e,e) + eT2(A,e,e), 



T,{A,e,e),Ai,o, ^(1 - ei^,),e), |(1 + £u;,)) 



(^A2{t, Ti{A, e, e), Ai,o, |(1 - SLO,),s) 




ti-- ('^3(t,Ti(yl,e,e) + er2(A,6»,e), 
f ^2,0, I (1 + eu;.), e), 63(4, Ti (A, 6, e) + 

'■® +ef2(Ae,£),A2,o,|(l + £a;,),£)) , 

V"^"---^ where ^2,0 = A2 (Ti (A, 0, e) + 

\ \ +£f2{A,L0,e),Tl{A,U),£), 

7r \ 

.X ^1,0, -(1 -£'^e),£ 

— > ^ ^ 

Acos ( Y^—e] 




/ \--4-r^smf-i-^6l') 



i (Ai(t,0,A,6i,£),6ii(t,0,A,e,e)) 



where Ai,o = Ai (ri(A, 0, e), 0, A, 6*, e) '^(i) (^iC^^iC^. " = e)= 0, A, (9, e), |(1 - £0;^)) 

Fig. 2: Illustration of sewing of the solution of system l l2.8t -l lZ8t by solutions of each its single component J2.8t . l |2.9t and l l2.10t . 



Summarizing, we can define tlie solution t ^ (^(t, ^, ^, e), ^(t, ^, ^, e)) of system (I2.8I )- (I2.10I ) as follows 
(see Fig. Ill) 



A{t,A,e,e) 

e{t,A,e,e) 



if te[0,Ti{A,e,e)], 



A2 [t, Ti{A, e, e), A{Ti{A, 9, e),A, 9, e), -(1 - eo;,), e 
92 (t, Ti{A, 9, e),A{Ti{A, 9, e), A, 9, e), ^(1 - euJe),e 

if te[mA^9,e),Ti{AJ,e)+ef2{A,9,j)], 
/As{t, Ti{A, 9, e) + eT2iA, 9, e), ^(Ti(A, 9, e) + eT2{A, 9, e), A, 9, e)\ 



-(1 + ew£),e 

^3(t, ri(A, e) + ef2{A, 9, e),A{Ti{A, 9, e) + eT2{A, 9, e),A, 9, e) 



TT 



\ 



if t e[Ti{A, 9, e)+eT2{A, 9, e),7r]. 

Step 2. At this step we show that fixed points of the Poincare map (12.111 ) can be studied by means of the 
function P introduced in the formulation of the theorem. To this end we decompose as 

'A^ 



P,{a, 0) = + e{P,,i{A, 9) + P,,2{A, 9) + Pe,z{A, 9)), 



where 



Ti(A,0,e) 



PeM,9) 



Gi{T,A{T,A,9,e),9{T,A,9,e),e)dT, 







Ti{Afi,e)+eT2{A,e,s) 



I 

'A,e,, 

TT 

I 



G2{T,AiT,A,9,£),9{T,A,9,e),£)dT, 



Ti(A,e,e) 

TT 



G3{T,AiT,A,9,e),9{T,A,9,£),e)dT. 

Ti{A,d,e)+eT2{A,e,e) 

Since sin, cos and g are bounded on any bounded set then from system (I2.8I) - (I2. 101) we have that 



A{t,A,9,e] 
9{t,A,9,e) 



A 
t + t 



as £ ^ 
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uniformly with respect to t G [0, vr], |j4 — Ao| < \6 — 6o\ < 6. This gives us immediately that 

_ Ti(A,6»,0) 

/ Gi{T,A,T + e,0)dT, 


P,,3(A^)^ / G3iT,A,T + e,0)dT, as e^O, 

Ti{A,e,0) 



(2.12) 



uniformly with respect to |A — Ao| < 6, \9 — 9o\ < 6. Since we proved that Ti and T2 are continuously 
differentiable, then (12.121 ) implies that 

(PeJ{A,e)^{Po,iy{A,e), 

{P,J{A,6)^{Po,3nA,9), as e ^ 0, 

uniformly with respect to |A — ^o| < ^, W — 0o\ < ^■ 

Let us now study the behavior of ^^^2 and {Pe,2)' as e — > 0. We have 



/ 1 



Ti{Afi,e)+eT2{a,e,e) 



-(1 - euje 



Ti{A,e,e) 



sm 



e \euje 
1 

A{T,A,d,e) 



^ 'e{T,A,e,e)-^{l-euJe)) + ''^ ^ 



1 



(1 - euJs)^^e COS {0{t, a, e, e) 



TT 



V 



{l-euoe)] + 



TT 



g r, £A{T,A,e,e) 



1 



1 - eujr. 



1 



TT 



COS ie{T,A,e,e)--{l-eu:e)]+-] , 



vr 



-A{T,A,e,e) 



1 - 



1 



TT 



sin 6'(r, A, 9, e) (1 - eus) + - , e Ur 



vr 



Scaling the time in the integral as r = Ti{A, 9, e) + eio^s we get 



T2(A,0.e)V, 



Pe,2iA,9) = -u;eil-euJe 



sin ( — (0(ri(A,6',e) +eswe, A,6l,e)- 



vr 



-(1-60;,) +- 
1 ^ 



vr 



. (1 — euj£)uj£- 

A{Ti{A,9,£) +eiOeS,A,uj,e) 

•cos ( — {9{Ti{A,9,e) +euJeS,A,9,e)- 



euJe: 



\ 



IT 



g Ti{A,9,e) +euj£S,eA{Ti{A,9,e) +euj£S,A,9,e)- — cos {9{Ti{A,9,e)+ 



1 — SLVp 



1 



IT 



+£0JsS,A,9,£) - -{l-eus)) ,-A{Ti{A,9,£)+eujsS,A,9,£) 



1 



1 — etOr 



1 



TT 



sin (9{Ti{A,9,£)+£uJeS,A,9,e) - -{1-euje)) +- ,e Us 



vr 



Put 



Since 



K,{A, 9) = - {9{Ti{A, 9, e) + ecoes, A, 9, e) - 9{n{A, 9, e),A, 9, e)) 



- {9{Ti{A,9,e)+euJeS,A,9,e) - ^{1 - euJe)] 
= Ke{A,9) ^ (9)[^){Ti{A,9,Q),A,9,Q)u:os = UQS, as e ^ 0, 



then 



vr 

sin I s H 



vr 



vr 



2 ^ 5 TT - ^) 0, ^ sin s + - , ds, as e ^ 
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uniformly with respect to \A — Aq\ < 6,\9 — 9q\ < 6. Since (Ke)'(j4,0) converges as e ^ uniformly in 
|A - ^ol < <^ and - 6'o| < then [KsYiA, 6) {Kq)'{A, 9) as e ^ 0. Therefore, 

(Pe,2y{A,e)^(Po,2yiA,e) as £^0 

uniformly with respect to |A — Ao| < 6,\9 — 9q\ < 6. 
Summarizing, we proved, that 



1 



Pe{A,t 



A 



PiA,9) + 



HPe,i - PoAA, 9) + P,,2(A, 9) - Po,2iA, 9) + Pe,z{A^ 9) - Po,3{A, 9)). 
Therefore, for any e > sufficiently small, the function 

{A,9)^P,{A,9) 

has a unique zero {As,9s) such that {As,9g) (^O)^o) as e ^ and the real parts of eigenvalues of 
{PcYiAg, 9g) — / are negative. This is equivalent to say that the eigenvalues of the matrix ((P£)'(^£, 9^))"^ be- 
long to the interval [0, 1). This implies (see 161 Lemma 9.2]) that 
t i—f {A{t, As,9^,e),9{t, A^,9s, e) — t) is an asymptotically stable vr-periodic solution of (I2.8I )- (I2.10I ). To 
see the latter one will probably wish to make the change of variables = ^(t) — t in (I2.8I )- (I2. 101 ). Since the 
change of variables (I2.5I )- (I2.7I ) is vr-periodic, than given by these formulas corresponding solution (xg, x^) of 
(I2.3I )- (I2.4I ) is also vr-periodic and asymptotically stable. Uniqueness of Xe follows from uniqueness of {A^, 9e) 
and the fact that the change of variables (I2.5I )- (I2.7I ) is one-to-one. 



To finish the proof it remains to observe that for any t G 

1 



vr 



^ A{t,Ae,9e,e) cos 



1 — euj, 



-9{t,A,,9e.e) 



\ 



-A(t,Ae,9e,£) 



1 - eijJr. 



■ sm 



^ 9{t,A„9„e, 



1 — £UJ. 



we have 



^0 cos(t + 6io) 
-Aq sin(t + 9q) 



as e 



and that for any t G [ — — 9q^tt 



vr 



we have 



V 



A{t,Ae,9e,£) COS 

-A{t,Ae,9e,e) sin 



1 — eojp 



1 - £L0^_ 



vr 



9{t,Ae,9e,£)- -{l+ELOe 



vr 



9{t,Ae,9e,e) - -{l+£UJe)] + - + 



vr 

+ -+vr 
vr 



as e 



0, 



Aq cos(t + 6^0 + vr) 
-Ao sin(t + 6*0 + vr) 

that is Xe converges to the solution of (11.21 ) with the initial condition ^^(O) = {Aq cos —Aq sin ^o) as e ^ 
pointwise on [0, vr]\{to}. The proof is complete. 



□ 



3 An application 

In this section we apply the result of section 2 to an impact oscillator shown in Fig. 3. 
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Fig. 3: Model of a preloaded ball (body of mass m = 1) bouncing 
against a nearly elastic surface S which is represented by a spring 

of stiffness — — - with e > small. The rest coordinate of the mass 
is assumed to coincide with the origin of the x-axis on the one hand 
and with the surface S at rest on the other hand. We also suppose 
that the viscous friction coefficient equals eci outside the contact 
with S and that it takes the value C2 + eci during the contact. 



//////////// 



A body of mass m = 1 is bouncing against a nearly elastic surface S (large stiffness l/e^cj^). Assuming in 
addition that the body is subjected to Rayleigh excitation, viscous friction and forcing, the equation of motions 
can be written as follows 

X + X = —eax — ecix + efiix{l — x'^) + sin t, if x ^ 0, 

1 (3.13) 

X -\ — 7^ — 7tX = —(c2 + £Ci)x + (fi2 + £tJ'i)x(l — x'^) + sint, if x<0. 

System (13.131 ) is of the form (11.11 ) with limiting system. There is no obvious reason why the Rayleigh excitation 
should be (9(1) during impact, but as this does not complicate the analysis, we admit this possibility. The 
averaging function P takes now the longer form 

-/iiAsin(T + 9){1 - sin2(r + 6*)) + 7sinr - 26j^cos(t + ff))dT - 
- (l cos(r + , + y (-«^ cos(r + ^ + .) + c,A sm(r + ^ + vr) + 

-^ii^sin(T + 6* + 7r)(l - £■ sin^(r + 6* + vr)) +7sinT - 2wAcos(t + 6* + 7r))(iT - 



-fO cos e - lA{ci + C2Cj) + f + 112^) A (1 - lA^) 
-^7cos0 - ^76'sin6' + f (a + 2uj) 



To formulate our result we need some preliminary notations. First, we introduce the function M : M — > M as 
follows 

^ / ,cos^ + 0sin^ 

M{e) = -OcosO + {ci + C2u:) + 

a + 2oj 



-{m + H2uj) 



cos 9 + 9 sin 9 87 / cos 9 + 9 sin 



a + 2io vr^ V a + 2uj 




Proposition 3.1 Let M {9o) = Ofor some 9q € (0, ^) and Aq is defined as 

7 cos 00 + 7^0 sin 6*0 
^0 = FT r7r\ • (3-14) 

If 

M'{9o) > (3.15) 

and 

-(ci + C2u;) + {m + ^i2u;)(l - 3Ao^) < 
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then, for any e > sufficiently small, equation ^3.13^ has exactly one n-periodic solution such that 

(rEe(O), i;e(0)) {Aq cos —^0 sin 6*0) as e ^ 
The solution Xe is asymptotically stable. 
The proof of proposition (13.11 ) relies on the following lemma. 



Lemma 3.2 Consider 2x2 real matrix D. If Sp D < and det ||D|| > then the eigenvalues of D have 
negative real parts. 

The statement of the lemma follows from the direct computation of the eigenvalues of D according to the stan- 
dard formula for roots of quadratic equations. 

Proof of proposition 13.11 Direct computation shows that {Aq, 6q) is a zero of P. To prove the proposition it 
remains to show that 

1) SpP'{Ao,eo) <0. 

2) det||P'(^o,^o)|| > 0. 

But these two relations follow from the formulas 

1) Sp P'{Ao, 60) = -7r(ci + C2w) + 7r(/ii + /U2w)(l - SA^). 

2) , ' . det\\P'{Ao,9o)\\=M'i9o), 
which ai^e straightforward. 

□ 

Our next proposition 13.31 shows that proposition 13. II is not vacuous, namely we give sufficient conditions 
ensuring that (13.151) is satisfied. Before proceeding to the formulation of proposition 13.31 we need to introduce 
some notations and properties. First we observe, that the 

Ocosd 

where 

KjO) = { ^1±^ - ^^1±^ fl - ^ / cosg + gsing y\\ A 
^ ' \ a + 2uj a + 2uj \ -k^ \ a + 2uj \9 



Observe, that there exists 9^ G (O, ^) such that 



vr 



K'{9) + for all g - j . (3.16) 

In fact, if K'{9n) = for some sequence ^„ | | as n ^ 00, then 

M'{9n) = {cos9n — 9nsm9n){—l + K{9n)) — > oo as n ^ oo, 
that contradicts the boundedness of the derivative of M. 

Proposition 3.3 Let 9^ £ (O, be such a number that Ii3.16\l holds true. Assume that 

K{9^) < 1 

and denote by 9q G (6**, f ) the unique point satisfying 

-1 + K{9q) = {). (3.17) 
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Then M{6q) = and M'{6q) < 0. Consequently, for any e > sufficiently small, equation ([3.13\l has exactly 
one ir-periodic solution Xs such that 

{xi;{0), X!r{0)) — > {Aq cos Oq, —Aq sin^o) as e ^ 0, 

where 

^ _ 7 cos 00 + 76*0 sin Op 
°~ f(a + 2u;) 

The solution is asymptotically stable. 

Note that ^ (^*) f ) satisfying (13.171 ) always exists and is unique since —1 + K{9^) < 0, K'{9) / 0, for 9 G 
(^*,f) ,and 

limK(6') = +00. (3.18) 

Proof of proposition I3.3[ First, observe that M{9q) = 9q cos 9o{—1 + K{9q)) = 0. Second, we have 

M'{9) = {cos9 -9sm9)L{9) + 9cos9K'{9) 

and so M'{9o) = K'{9o). But properties (13.161 ) and (13.181 ) imply that K'{9q) > and so the proof is complete. 

□ 

Let us finally formulate our result in the simpler setting when the Rayleigh excitation is switched off, that 
is /ii = /i2 = 0. We have that 



a + 2io \9 

in particular, there exists an unique 9^ G (O, ^) such that sign K'{9^) = sign [9 — cos 9) = 0. 
Proposition 3.4 Let 9^ G (O) f ) be the unique point such that 0* — cos 9^ = 0. Assume that 

C1+C2UJ f 1 \ / 1 

■TT + tg 9* ] < 1. 



a + 2uj \9 

Denote by 9q G (^9^:, ^) the unique point such that 



a + 2uj \9q J 
Then the conclusion of proposition \3. 3\ holds true. 

4 Discussion 

• It is remarkable that the analysis of system (11.11 ) can be handled by the introduction of a Poincare map 
and the use of the implicit function theorem although the limiting system for e ^ is non-smooth. 

• At the same time we have formulated an unusual type of singular perturbation problem. Putting e = 0, 
we have non-smooth impact, for e > we have fast motion in a neighborhood of the subset x = 0. For 
a:: > slow motion takes place but this is not described by standard slow manifold theory, see ifTTI . Still, 
the dynamics for x > can be considered as taking place in an explicitly formulated slow manifold. 
On the other hand, the solutions for x < have as slow manifold the boundary x = 0. This does not 
satisfy the necessary hyperbolicity condition, but the solutions for x > are forced to the manifold x = 
and, after passing by a fast transition through the domain x < they ai^e forced again to leave x = 0. 
We note also that sliding along the slow manifold, as happens for instance in dry friction problems, is not 
possible. This simplifies the bifurcational behavior. 

• Regarding the averaging result obtained in this paper, we draw attention to the papers lITOl - 191 and 
further references there. In ifTOl the framework of differential inclusions is used, in [9J explicit estimates 
of the vector field and the solutions are given in the case of impulsive forces. Our approach economically 
avoids the estimate of general solution behavior as we aim at the more restricted result of obtaining 
periodic solutions. 
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